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Abstract 


We suggest a convenient method based on the Fibonacci polynomials and 
the collocation points for solving approximately the Abel’s integral equation 
of second kind. Initially, the solution is supposed in the form of the Fibonacci 
polynomials truncated series with the unknown coefficients. Then, by placing 
this series into the main problem and collocating the resulting equation at 
some points, a system of algebraic equations is obtained. After solving it, 
the unknown coefficients and so the solution of main problem are determined. 
The error analysis is discussed elaborately. Also, the reliability of the method 
is quantified through numerical examples. 
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1 Introduction 


It is well known that the Abel’s integral equation is an applied Volterra 
integral equation with a singular kernel. Some applications of it appear in 
the plasma spectroscopy [3], thermoelectricity [8], and spectrographic data 
[9]. 

Solving the integral equations with singularity feature is often contestable 
because of the complexity and the heavy computations. Thus, various meth- 
ods especially in recent years have been proposed for solving the Abel’s inte- 
gral equation. Some of them include Bernstein polynomials [2], Taylor expan- 
sion [10], quadrature and mechanical quadrature methods [11, 12], Laplace 
transform [16], Block-pulse functions [24], modified Homotopy and Adomian 
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methods [25], variational iteration method [26], shifted Legendre collocation 
method [27], Laguerre polynomials [28], and Legendre wavelets [29]. 


Recently, the nonorthogonal polynomials have been used widely to solve 
different problems; sce [5, 4, 14, 23]. One of these polynomials is the Fibonacci 
polynomials. They have been derived successfully to find the approximate 
solutions of the fractional diffusion equation [6], the boundary value prob- 
lems [15], the Fredholm-Volterra integral equations [19, 21], the systems of 
high-order linear Volterra integro-differential equations with variable coeffi- 
cients [22], the nonlinear stochastic It6-Volterra integral equations [20], the 
singularly perturbed differential-difference equations [17], and the systems of 
linear Fredholm integro-differential equations [18]. Moreover, some advan- 
tages of the Fibonacci polynomials rather than orthogonal polynomials have 
been described in [1]. 


Throughout this manuscript, the generalized Abel’s integral equation of 
second kind is considered in the following form: 


* J(u(t) 
u(x) = g(x) +k fo at, (1 
o (et) 
so thatO <a<1,k€ER,0<2,t < b for b > 0, and g and J are real 
continuous functions. Note that if a = 4 and J(u(t)) = u(t), then (1) is 
reduced to the conventional Abel’s integral equation of second kind. 


To the best of our knowledge, the existing works on the Fibonacci poly- 
nomials in the literature have mainly concentrated on the problems with 
nonsingular kernels. In this paper, we try to assess the application of these 
polynomials for solving an important family of the integral equations with the 
singular kernel. In fact, the aim of this research is to present an approximate 
solution of the problem (1) in the form of 


N 
u(x) % uy(2) = >> Gn Fn(2), (2) 
n=1 


where F,, n = 1,...,N denote to Fibonacci polynomials; a,, n = 1,...,N 
are the unknown coefficients; and N is chosen any positive integer. Here, the 
Fibonacci polynomials are defined by [7, 13] 


(°#] 


F,(2) = > Cae ee n=1,...,N. (3) 


Furthermore, the Fibonacci polynomials can also be identified by the recur- 
rence formula below 
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Fi, (2) = cF,_-1(x) + Fy_o(x), n =3,4,.... 


To achieve the aforesaid aim,(2) is substituted into (1) and by collocating the 
obtained identity at the collocation points x; = wal —1),j =1,...,N, 
a system of algebraic equations consisting of N unknowns is gained whose 
solution is a,,...,any. Therefore, after solving this system the solution of 
corresponding problem can be found by applying (2). 


In the next section, the Fibonacci collocation method for the problem (1) 
is formulated by utilizing the matrix forms. The error of suggested method 
for (1) is discussed in Section 3. Section 4 provides the numerical results of 
applying the Fibonacci collocation method by four illustrative examples. At 
the end, a conclusion is given in Section 5. 


2 Method of solution 


At the beginning, let us exhibit (3) in the matrix form as follows 


where 


(¢) oO oO 90 0 
0 () oO 0 0 
? e) : Qo 
D= : ‘ 
(2) 0 GA) 0 0 
0 (wi) 0 (xfs) MOD} 


if N is even, and 
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(3) 0 0 O: gee D 
0 (j) 0 0 cee OD 
G) 0 (-) 0 0 
p-| ® O © M0 
N-1 oe 
0 (ws) 0 (w2s) 0 
N-1 N+1 _ 
(is) 0 (a#s) co 


if N is odd. 
We divide the structure of the proposed method into linear case and 
nonlinear one. 


e Linear case: 


In this case, J(u(t)) = u(t). For the sake of simplicity, (1) may be 
shown in the form of 


u(x) = g(a) + kI(x), (5) 


in which 
I(x) = | : wt (6) 


For a given N, the desired solution of (1) is supposed in the form of the 
truncated Fibonacci series (2). So, u(x) can be written in the matrix 
form 

u(z) =F(a)A; A= |[a, a... ay], 


or with the help of (4), it follows that 


u(x) = X(x)D7A. (7) 


By substituting (7) into (6), it is obvious 


9= fara (f eat) = 


To acquire a matrix relation for I(x), an explicit formula must be de- 
termined for the integral 


Tio(& = —— Ht, FSO Aji ND 
a(2) | Ete 


For this purpose, changing the variables by t = xz enables one to write 
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dt=adz, 0<z<1, 


and hence 


1 
Tia(£) = aad) (1 — z)-%2*dz = B(i+1,1-a)2*t*"%, (9) 


0 


or equivalently 


eel 8) aap 
Kiel®) = “TG a +0) os 


where I'(-) and 4(-,-) refer to the Gamma and Beta functions; see [5]. 
Lastly, using (8) and (9) results in 

I(x) =I,(x)D7 A, (10) 
where 


I(x) = [B(1, 1 — a)a!-* B(2,1—a)a?-@ ... B(N, 1 — a)2X-9]. 


In order to implement the numerical procedure, by substituting (7) and 
(10) into (5) it is clear that 


X(x)D* A = kI,(x)D7 A + g(x). (11) 


Now, collocating (11) at the points 


b 2b 
{e1.tata,.. tw} = {0 op ee Ph (12) 
the main matrix equation is produced as follows: 
{XD* —kI,D’}A=G, (13) 
in which 
X(21) I, (1) g(x1) 
_ X(x2)} = I,(%2) g(x2) 
X= ’ I, = i ) G= 
X(zy) I.,(xn) g(zN) 


It is notable that A and G are column vectors with N entries and the 
dimension of others in (13) is N x N. 


Briefly, (13) may be shown as 
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WA =G, (14) 
so that 
W = XD? —&I,D’. 


Unequivocally, (14) is a linear system of algebraic equations whose so- 
lution is the coefficients a1, a2,...,ay. Finally, the solution of (1) can 
be uncovered by employing (7). 

Nonlinear case: 


Based on (7), we have 
u(a) = X(x)B, (15) 


where B = D7 A. If m € N and £, = (N — 1)m +1, then there exists 
the vector P,, = [pi pz... Pém| such that 


u™ (x) = Xm(x)Pm, (16) 
in which 


Xp (x) = [1 a ge ; 


According to the multi-binomial theorem, the entries of P,, are ex- 
pressed in terms of a1, @2,...,a@y. In the other word, for 1 <i, < Em, 


Pim = Pin, (Q1,42,-..,@n), where y;,, is a real multivariable function 
for every m and i. 


Suppose that J applied in (1) is an M € N times differentiable function 
at x = 0. Then, 


a m!) 
Therefore, 
M 
J) (0) sg 
T(u(t)) = S° u(t) (17) 
m=0 : 


Looking at (16) and (17), we conclude 


Mam) 


From (1), (15) and (18), we obtain 
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J(™) (0) 
X(x)B g(x) +k S$ Rp (a)Pm- (19) 


It must be mentioned that 


Rm(x) = [rili<i<en = [ at dt, 


x —t)% 
is a row vector and by using (9), one can imply that 
ri = B(i+1,1—a)a't*?, t=1,2,...,€m- 


If we replace ~ by = in (19) and collocating the obtained identity at 
the points defined by (12), a nonlinear system of algebraic equations is 
produced which may be solved by any available iterative scheme such 
as the Newton—Raphson method. 


3 A discussion on the error 


Throughout this section, the error of the method proposed in the previous 
section is analyzed. Consider L°[0, b] is the L-infinity space defined on (0, }] 
with the norm 


[P| = sup |h(x)|, 
O<2<b 
where fh is a real function. In addition, for F;(xz), 1 < 7 < N, suppose that 
F = span{F (x), Fo(x),..., Fn(x)} Cc L® (0, b]. 


Since F is a finite-dimensional vector space, h(x) has the best approximation 
such as h,(x) € F, namely 


for all r(x) € F: ||h(x) — hs(£)loo < |[A(x) — r(x) |loo- 


Due to h(x) € F, there exist the unique coefficients c, such that 
N 
h(a) © hs(2) = S> cp F(x) = F(x)C7, 
k=1 


where C = [er C2... cn| ; 


Lemma 1. Let h : [20,6] > R is N times continuously differentiable for 
0 <2 < b. If hs(x) is the best approximation to h, then the error bound is 
declared as follows 
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M(b — 20) 


(20) 
in which M = supg<x<p [A (x)|. 


Proof. Suppose that h(a) is an arbitrary approximation of h(x). We choose 
this approximation in the form of Taylor series of h(x). Clearly, 


as x — Xo 2 
h(x) = h(a) + h' (xo) (x Xo) + h" (ao) $ my ) 
~4 (2 —20)*— 
fee pA (80) TR DT 


Therefore, there exists 7 € (xo, b) such that 
j, _ |p), = 20) 
|h(x) — A(x)| = ph" (n) —=3—— 


Since h(x) is the best approximation of h, 


l|n(x) — hs(2)lloo S (wx) — A(x) loo = eeu. aa) = h(x)| 


N 
= (N) (x — £0) 
= sup [A\”’(n) ni : 


O<a<b 


Putting M = supo<r<p [A (x)|, (20) is realized. 


The main result of this section is announced in the following theorem. 


Theorem 1. Assume that u(x) and uy (x) for a given N > 1 are the exact 
and the Fibonacci solutions of (1), respectively, on [0, 6]. If J satisfies in the 
Lipschitz condition with the constant 7 and also 1 — a — |k|nb'~* > 0, then 
one can deduce 


M'(b — xo)N 


I|u(x) — un ()|loo S (21) 


1—a—|k|nbt~% 
l-a 


where \ = and M’ = supg<e<s |g) (2)|. 


Proof. Suppose that g(x) is expanded in terms of the Fibonacci polynomi- 
als; then the obtained solution is an approximated polynomial, uy (x). We 
attempt to seek an upper bound for the associated error between the exact 
solution u(x) and the approximated solution uy (x) for (1). Consequently, 


u(x) — un (x) = g(x) — gn (a) 4 kf Ha) ~ T(un(t)) 5 


Applying the triangular inequality, it is concluded 
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* J(u(t)) — 
jue) — un (@)] < fe) — aay + | f° AP 
Besides, one can write 


[ J(u(t)) = ae < \| J(u(a)) _ J(un(2))|loo A dt 


(x — t)@ (x — t) 


= FT ll4(u@)) - Jun @))Iloo: 
Thus, for x € [0,6], we have 
* F(u(t)) — Tun) yl < ne 
fe O'all <P tule) — unlaliw (23) 


On the other hand, if M’ = sup |g) (x) for x € [0,b], based on (20), then 


M'(b—a0)% 
lo(2) — gn(a)| < Ot (24) 
Ultimately according to the relations (22)-(24) and the mentioned assump- 
tion, (21) is satisfied. 


4 Numerical examples 


This section is devoted to four examples for the description of the method’s 
applicability. The corresponding computations are performed by MATLAB 
R2015a software on a 64-bit PC with 2.20 GHz processor and 8 GB memory. 


Example 1. [16] Consider the Abel’s integral equation (1) in the form 
of 


The exact solution of (25) is u(x) = x?. 


Suppose that N = 3. By this selection for N, the collocation points 


become {0, $ 1}. After some calculations, we have 
1 0 0 0 
W = |985/408 1155/1189 1027/360] , G= | (2/2/15) + (1/4) 
3 7/3 76/15 31/15 


By solving the algebraic system WA = G, one obtains 
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A = [-1 —5.67 x 10726 1)”. 


Accordingly, u3(x) = (—5.67 x 10~'®)x + x?. This solution concludes that 
\|u(x) — u3(x) loo = 5.67 x 1071°. It is worth to mention that in [16], the 
maximum absolute error of the Laplace transform method is 5 x 107” at the 
level n = 25, which requires large computational effort. 


Example 2. [25] Consider the Abel’s integral equation (1) as follows: 
" t 
u(x) = 2/x — | tO) gy. nui, (26) 
o (a—t)?2 


The exact solution of (26) is u(x) = 1—e”*erfc(./7a), in which the comple- 
mentary error function erfc is defined as 


erfc(x) = a i. ew dr. 


By deriving the presented method for N = 5 and N = 10, the approximate 
solutions are 


us(x) = —5.257011824 + 12.97393123 — 11.381689x? + 4.37816712 
+(2.8421709 x 10-14), 

uio(x) = 318.98646x° — 1582.07372% + 3371.41952" — 4038.71492° 
+2986.89022° — 1409.7584a4 + 424.24302xr3 — 79.328853x7 + 9.0564174z 
+(9.3132257 x 1071°). 


Table 1 and Figure 1 show the numerical results and the absolute errors of 
the presented method, respectively. It is seen the results of applying the 
proposed method are close to the exact solution. One should note that the 
corresponding values of the exact solution have been calculated with the aid 
of erfc(-) command of MATLAB software. 


Table 1: Numerical results of Example 2 
Presented method 


t Wos V=10 Exact solution 
0.0 0.00000 0.00000 0.00000 
0.2 0.51575 0.51305 0.50835 
0.4 0.62595 0.60539 0.60335 
0.6 0.65055 0.65762 0.65632 
0.8 0.70763 0.69299 0.69184 
1.0 0.71340 0.71974 0.71794 


Example 3. [24, 29] Consider the following Abel’s integral equation 


The Fibonacci polynomials solution for Abel’s integral ... 73 


Figure 1: The illustration of absolute error for Example 2 


_t 1 Lat l-«z 1 f* u(t) 
al A ra nl (=) al Te (27) 


The exact solution of (27) is u(a) = Jat 


The solutions of the Fibonacci collocation method for N = 4 and N = 8 
are identified as follows: 


ua(x) = —0.0852864x°* + 0.277508x? — 0.4851932 + 1.0, 


ug(x) = —0.01209192" + 0.06155552° — 0.1473522x° + 0.23416724 
—0.303006x? + 0.373771a? — 0.499937a + 1.0. 


Figure 2 is connected with the absolute error of the approximate solutions 
in the logarithmic scale. In Table 2, the numerical results of the suggested 
method have been compared with those of Block-Pulse functions (BPFs) 
method for 32 and 64 basis functions [24] and the Legendre wavelets (LWs) 
method for m = 10 ones [29]. It seems the presented method provides more 
accurate results than BPFs and LWs methods by smaller number of basis 
functions. 


Example 4. [30] As the final example, consider that 


u(x) = [ 7 -u(t)dt + (2), (28) 


3 _ 4096.17 3 


in which g(x) = x x2. The exact solution of (28) is u(x) = x”. 
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Table 2: Numerical results of Example 3 

BPFs [24] LWs [29] Ours 
k = 32 k= 64 m = 10 N=4 N=8 
0.0 0.999123 0.999993 0.999432 1.000000 1.000000 1.000000 
0.2 0.912305 0.912873 0.912320 0.913379 0.912871 0.912871 
0.4 0.845156 0.845154 0.845321 0.844866 0.845154 0.845154 
0.6 0.790527 0.790562 0.790539 0.790365 0.790569 0.790569 
0.8 0.745361 0.745316 0.745342 0.745784 0.745356 0.745356 
1.0 0.707120 0.707103 0.707163 0.707029 0.707107 0.707107 


Exact 


Figure 2: The illustration of absolute error for Example 3 


For solving (28) by the presented method, similar to the previous examples 
assume that the approximate solution is in the form of (7). Subsequently, 


u?(«) = X(x)D7 AX(z)D* A = A?DX" (x) X(x)D7A. (29) 
Placing (29) into (28), it is obvious that 


X(t) 
Va-t 


X(c)D?A = rATD O a) DA + 9(2). 


Note that X(t) = ¢X7(t)X(t) = [ajli<ij<n is an N x N matrix so that 
cjj = t'tI-1. With the aid of (9), we can write 
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Q(z) = [gijhi<ij<n = 2 Uf Joa 


in which 
: ‘ 1 i+j+2 
gj =B\ its, )ar?. 
Hence, the following main matrix equation is found 
X(xz)D7 A = A7DQ(z)D7A + g(z). (30) 


Collocating (30) at the points defined as (12) generates a nonlinear algebraic 
system consisting of N equations with N unknowns. We solve this system 
by use of MATLAB fsolve command with the initial guess 


Aq = (00... 0]. 
N 


Clearly, when A is detected the solution of (28) can be assigned by (7). 


Utilizing the described method for N = 5 and N = 8, the solutions are 
acquired as follows: 
us(x) = —0.000027344159x* + 1.0000417x* — 0.000019310296x? 
+0.0000026504592z + (3.3881318 x 10~?'), 
ug(x”) = —0.00019083629z7 + 0.00055866874x° — 0.00064855817x° 
+0.00038192279x* + 0.99987881a° + 0.00001972285627 
—0.0000012761595a + (1.0842022 x 10719). 


Figure 3 is related to the absolute error of the obtained solutions in the 
logarithmic scale. 

Table 3 compares the absolute errors results of the proposed scheme for 
N basis functions with the second kind Chebyshev wavelet method [30] for 
m! = 2*-1M ones at some points of [0,1]. It is observed that the presented 
method produces more accurate solutions. 


5 Conclusion 


In this research, a new method using the Fibonacci polynomials, collocation 
points, and matrix operations was established for solving the Abel’s integral 
equation of second kind. The useful properties of the proposed method with 
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Figure 3: The illustration of absolute error for Example 4 


the help of collocation technique converted the problem to a solvable system 
of algebraic equations. As it was seen, to gain the best approximations, the 
value of N must be selected large enough. 

One of the method’s advantages is the matrix forms of the calculations. 
This style of writing makes the computer programs more comfortable. The 
examples demonstrated the validity and the correctness of the presented 
method by some tables and figures. The numerical results and the compar- 
isons with other methods in the literature revealed that the presented method 
works well in terms of the suitability, impression, accuracy, and computer- 
orientation. 

The method described in this paper may be developed for other problems. 
Especially, the fractional integro-differential equations with a weakly singular 
kernel or the systems of Abel’s integral equations can be suitable problems 
for everyone who is interested in these areas. 
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